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AERONAUTIC  SYMBOLS 
1.  FUNDAMENTAL  AND  DERIVED  UNITS 


w 

g 

m 

I 

n 


S 

G 

b 

c 

A 

F 

2 

L 

D 

Oo 

Di 

Dp 

C 


Symbol 

Metric 

English 

Unit 

Abbrevia¬ 

tion 

Unit 

Abbrevia¬ 

tion 

Length _ 

1 

meter _ _ _ 

m 

foot  (or  mile) _ 

Time _ 

t 

second _ ; _ 

8 

second  (or  hour) _ _ 

Force _ 

F 

weight  of  1  kilogram _ 

kg 

weight  of  1  pound _ 

lb 

Power _ 

P 

horsepower  (metric) _ 

horsepower..  . 

fkilom'eters  per  hour _ 

kph 

miles  per  hour _ 

mph 

Speed - 

V 

(meters  per  second _ 

mps 

feet  per  second _ 

fps 

2.  GENERAL  SYMBOLS 


Weight=m^ 

Standard  acceleration  of  gravity=9.80665  m/s’ 
or  32.1740  ft/sec’ 

W  '  , 

Mass=-j  •  .  -  ,  ■  .  ■ 

Moment  of  inertia=mP.  (Indicate  axis  of 
radius  of  gyration  k  by  proper  subscript.) 
Coefficient  of  viscosity 


p  Kinematic  viscosity 

p  Density  (mass  per  unit  volume) 

Standard  density  of  dry  air,  0.12497  kg-m~'‘-s*  at  15' 
and  760  mm;  or  0.002378  Ib-ft"*  sec’ 

Specific  weight  of  “standard”  air,  1.2255  kg/m’ 
0.07651  Ib/cu  ft  '  i. 


3.  AERODYNAMIC  SYMBOLS 


Area  ", 

Area  of  wing 

Gap 

Span 

Chord 

Aspect  ratio,  ^ 
True  air  speed 


1 


Dynamic  pressure, 

Lift,  absolute  coefficient  Cl= 


Drag,  absolute  coefficient  Cd= 


L 

qS 

D 


qS 


Do 


Profile  drag,  absolute  coefficient 
Induced  drag,  absolute  coefficient 

*  qo 


Parasite  drag,  absolute  coefficient  Ccp= 


qS 


Cross-wind  force,  absolute  coefficient  Cc= 


G 

qS 


it- 

Q 

Q 

R 


.  a 
€ 


da 


Angle  of  setting  of  wings  (relative  to  thrust  li 
Angle  of  stabilizer  setting  (relative  to  thr 
line)  ,  . 

Resultant  moment 

Resultant  angular  velocity  . 

Reynolds  number,  p—  where  Z  is  a  linear  dim 

sion  (e.g.,  for  an  airfoil  of  1 .0  ft  chord,  100  m 
standard  pressure  at  16°  Cj  the  correspond 
Reynolds  number  is  935,400;  or  for  an  air 
of  1.0  m  chord,  100  mps,  the  correspond 
Reynolds  number  is  6,865,000) 

Angle  of  attack  : 

Anglo  of  downwash 
Angle  of  attack,  infinite  aspect  ratio 
Angle  of  attack,  induced 
Angle  of  attack,  absolute  (measured  from  ze 
lift  position) 

Flight-path  angle 
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IN  A  SUPERSONIC  MAIN  STREAM 

l$y  I.  li.  Garhick  and  S.  I,  Rubinow 


SUMMARY 

A  theoretical  study,  based  on  the  linearized  equations  of 
motion  Jor  small  disturbances,  is  made  oj  the  air  forces  on 
wings  of  general  plan  forms  moving  forward  at  a  constant 
supersonic  speed.  The  boundary  problem  is  set  up  for  both 
the  harmonically  oscillating  and  the  steady  conditions.  Two 
types  of  boundary  conditions  are  distinguished,  which  are 
designated  “purely  supersonic’  and  “mixed  supersonic.”  The 
purely  supersonic  case  involves  independence  of  action  of  the 
upper  and  lower  surfaces  of  the  airfoil  and  the  present  analysis 
is  mainly  concerned  with  this  case.  .1  discussion  is  first 
given  of  the  fundamental  or  elementary  solution  corresponding 
to  a  moving  source.  The  solutions  for  the  velocity  potential 
are  then  synthesized  by  means  of  integration  of  the  fundamental 
•solution  for  the  moving  source.  The  method  is  illustrated  by 
applications  to  a  number  of  examples  for  both  the  steady  and 
the  oscillating  cases  and  for  various  plan  forms,  including 
swept  wings  and  rectangular  and  triangular  plan  forms. 
The  special  result's  of  a  number  of  authors  are  shown  to  be. 
mcluded  in  the  a:mlysis.  .  ■  'V  ■  •  :  y 


This  paper  constitutes  a  theoretical  study  of  the  aero¬ 
dynamic  forces  on  an  oscillating  or  steady  wing  of  finite 
span  moving  forward  at  a  uniform  supersonic  speed.  The 
treatment  is  based  on  the  linearized  theory  obtained  by 
considering  only  small  disturbances  in  an  ideal  fluid.  The 
wing  is  therefore  considered  to  be  a  nearly  fiat  thin  surface 
at  a  small  angle  of  attack  and  the  flow  is  considered  non- 
viscous  and  free  of  strong  shocks.  The  theory  in  this  case  is 
equivalent  to  finding  certain  solutions  of  the  wave  equation  in 
three  dimensions  with  respect  to  a  moving  coordinate  system. 

For  the  case  of  steady  motion  there  e.xist  a  number  of 
interesting  solutions  and  methods.  Among  these  may  bo 
mentioned  the  Von  Kiirman  and  Moore  linearized  treat¬ 
ment  of  slender  bodies  of  revolution  (reference  f),  the  Prandtl 
acceleration-potential  method  employed  by  Schlichting 
(references  2  and  3),  the  Busemann  method  of  “linearized 
conical  flows”  (reference  4),  studies  of  Jones,  Puckett, 
Stewart,  Brown,  and  Gurevich  (references  5  to  9);  and  a 
method  of  Von  Kdrmdn  employing  Fourier  integral  solutions 


of  the  two-dimensional  wave  equation  and  described  by 
him  as  “acoustic  oscillator  method”  (Wright  Brothers 
Memorial  lecture,  Dec.  17,  1946). 

The  corresponding  unsteady  or  nonstationary  problem  for 
two-dimensional  flow  (infinite  aspect  ratio)  may  be  con¬ 
sidered  to  be  solved.  In  this  connection  there  may  be 
mentioned  the  work  of  Possio,  Von  Borbely,  Temple  and 
Jahn,  and  the  present  authors  (references  10  to  13).  Of 
interest  also  are  two  wartime  German  papers  by  Schwarz 
and  Honl  (references  14  and  15).  The  corresponding  steady 
plane  case  to  which  the  nonstationary  problem  may  be 
reduced  is  that  treated  by  xickeret. 

Results  for  the  nonstationary  or  oscillating  case  are  of 
great  interest  in  the  investigation  of  aircraft  instability. 
The  two-dimensional  results  have  been  applied  to  a  study  of 
flutter  at  supersonic  speeds  in  references  12  and  13.  Of 
more  direct  interest  for  this  application  are  the  three- 
dimensional  results,  especially  for  wings  of  swept  plan  form. 

The  method  used  in  the  present  study  is  to  build  up  solu¬ 
tions  of  the  equation  satisfied  by  the  velocity  potential  by 
superposition  of  the  fundamental  wave-potential  solution 
for  a  spherical  source.  These  solutions  are  also  made  to 
satisfy  certain  required  boundary  conditions  on  the  airfoil 
surface.  In  the  two-dimensional  supersonic  nonstationary 
case,  which  appears  herein  as  a  special  limiting  case,  it  can 
be  proved  tliat  the  i)rocedure  leads  to  a  solution  that  is  the 
unique  solution  of  the  given  boundary  problem.  (For  the 
problem  of  subsonic  flow  past  a  thin  wing,  reference  may  be 
made  to  the  general  treatment  and  method  of  Kiissner 
(reference  16)  which  also  involves  solutions  of  the  wave 
equation.) 

Some  (jualitative  features  of  the  nature  of  the  boundary 
problem  may  be  mentioned  here.  F urther  remarks  may  be 
found  in  reference  17  and  in  Von  Kiirmiin’s  Wright  Brothers 
Memorial  lecture.  In  the  case  of  subsonic  flow  past  an  air¬ 
foil  the  whole  field  is  influenced  by  the  body.  The  concept 
of  circulation  has  proved  to  bo  very  useful  and  the  Kutta 
condition  has  been  used  to  specify  the  circulation  by  requir¬ 
ing  smooth  flow  leaving  the  trailing  edge.  Thus,  a  deflected 
aileron  in  subsonic  flow  influences  the  flow  pattern  over  the 
whole  wing  even  more  importantly  than  over  the  aileron  itself. 
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In  the  case  of  supersonic  flow  the  influence  of  the  body  is 
limited  to  only  certain  parts  of  the  field  of  flow  and  generally 
the  wake  does  not  influence  the  upstream  flow.  The  bound¬ 
ary  problem  for  a  three-dimensional  surface  moving  at  a 
supersonic  speed  can  be  classified  into  two  types  referred  to 
herein  as  “purely  supersonic’’  and  “mixed  supersonic.” 
The  definition  of  these  terms  is  given  in  the  analysis  accord¬ 
ing  to  the  parts  of  the  field  influenced  by  the  airfoil,  the 
purely  supersonic  case  involving  independence  of  action  of 
the  top  and  bottom  surfaces  and  no  reflecting  surfac'es  in  the 
field.  Thus,  in  the  purely  supereonic  case,  a  deflection  of 
the  aileron  would  produce  oidy  a  local  efl'ect  near  the  aileron; 
in  the  mixed  supersonic  case,  it  may  have  a  decided  influence 
on  the  jiart  of  the  wing  adjacent  to  the  aileron  or  on  other 
parts  of  the  wing.  For  a  given  wing  both  types  of  problems 
may  bo  involved. 

The  treatment  used  for  the  purely  supersonic  cases,  involv¬ 
ing  source  and  sink  distributions  to  account  for  the  action 
of  the  body,  is  believed  to  be  exact  within  the  framework  of 
the  linearized  theory.  The  upper  and  lower  surfaces  of  the 
airfoil  are  regarded  as  acting  independently,  each  surface 
being  “unaware”  of  the  presence  of  the  other.  The  treat¬ 
ment  is  thus  analogous  to  that  of  sound  in  a  moving  medium 
generated  by  the  motion  of  pistons  imbedded  in  an  infinite 
plane.  This  flow  picture  is  obviously  incomplete  in  the 
mixed  case  and  more  complicated  distributions  (doublets) 
are  also  required.  For  some  purposes,  however,  the  simpler 
treatment  may  still  be  used  in  conjunction  with  appropriate 
correction  factors.  Also,  for  steady  flow  past  a  symmetrical 
airfoil  at  zero  lift,  the  simpler  treatment  can  be  employed 
for  study  of  the  wave  drag. 

The  object  of  the  present  paper  is  to  develop  the  expres¬ 
sion  for  the  velocity  potential  in  the  purely  supersonic  case, 
based  on  the  elementary  solution  for  the  sound  source 
moving  uniformly  at  a  supersonic  speed,  and  to  indicate 
its  application  by  a  number  of  special  examples. 

SYMBOLS 


<t> 

disturbance-velocity  potential 

rectangular  coordinates  for  fixed  system 

x,y,z 

rectangular  coordmates  attached  to  source  mov¬ 
ing  in  negative  x-direction;  also  represents 
field  point  being  influenced 

i'.n'A' 

rectangular  coordinates  used  to  represent  space 
coordinates  in  fixed  system 

i,n,i 

rectangular  coordmates  used  to  represent  space 
location  at  source  distribution  xl({,  tj,  f) 

t,  T,  t' 

time 

V 

velocity  of  main  stream 

c 

velocity  of  sound 

M 

Mach  number  (r/c) 

r 

distance  defined  by  equation  (8) 

Tl,  T2 

time  function  defined  in  equation  (7a) 

1 

g 

function  defining  airfoil  surface  (y=g(_x,  z,  t)) 

fo,  u,  i'-i, 

e 


V 

To 

p 

a 

a 


W(J,  c,  t) 


A 

h 

li 


limits  defined  in  equation  (10) 
variable  used  instead  of  f  defined  by 
preceding  equation  (15a) 
pressure 

reference  pressure 
density 

angle  of  attack 
time  derivative  of  a 
angular  frequency 

vertical  velocity  factored  in  equation  (I 
space  function  ir(r,3)  and  time  fund 
angle  of  sweep 
vertical  displacement 
time  derivative  of  h 


ANALYSIS 

WAVE  EQUATION  AND  SOURCE  SOLUTIONS 

In  the  linearized  theory  based  on  small  disturbai 
equation  satisfied  by  the  velocity  potential  for  the 
gation  of  sound  waves  of  small  amplitude  is  the  wave  C' 


The  fluid  rhedium  is  considered  at  rest  at  infinity. 

In  the  treatment  of  linear  partial  differential  equati' 
so-called  elementary  or  fundamental  solution  is  o 
importance  since  general  solutions  can  be  built  up 
tributions  of  elementary  solutions.  From  a  physical  ] 
view  the  elementary  solution  may  correspond  to  a  sou 
discussion  of  the  nature  of  elementary  solutions  for  hyi 
differential  equations  of  a  general  type  has  been  gi 
Hadamard  (reference  18),  who  makes  the  cardinal  sta 
that  “every  result  of  the  theory  can  be  and  has  to  be  d 
from  the  consideration  of  the  elementary  solution  on 

A  fundamental  solution  of  equation  (1)  from  which 
solutions  may  be  formed  is  that  of  a  source  of  soun 
in  the  medium 


where 

r' = ^i{x'-ir+{y'-v'y+{z'-i'y 


In  equation  (2)  the  fixed  source  is  located  at  the 
({',  f'),  the  strength  of  the  source  is  A({',  j 

and  the  minus  sign  indicates  that  the  spherical  wa\ 
ili verging  from  the  center  of  the  disturbance. 

Another  closely  related  solution  of  equation  (1)  is  ; 
a  fixed  point  source  for  which  the  spherical  waves  ai 
verging  onto  the  source 


The  wave  potential  in  equation  (2)  is  often  desi; 
“retarded”  and  that  in  equation  (8),  “advanced.” 
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It  is  intoiulod  to  considi'i-  tliin  lifting  siirfacos  of  small 
cui  vaturc  which  arc  moving  forward  at  a  constant  supersonic 
velocity  r  and  which  may  he  performing  small  oscillations 
normal  to  the  direction  of  r.  The  direction  of  v  will  be  that 
of  the  negative  r-axis  and  the  surface  will  be  replaced  by  a 
distribution  of  moving  sources  in  the  is-plane  (lig.  1). 


y 


Figure  l.— 'I'he  coordinate  system  used.  Thin  lifting  surface  in  ir*plano  moving  at  a 
constant  supersonic  velocity  r  in  the  negative  x-dircction. 


Consider  a  source  moving  in  the  negative  x-diroction  with 
uniform  velocity  v  and  a  rectangular  coordinate  system 
attached  to  the  moving  source.  If  the  new  coordinates  are 
designated  by  x,  y,  z,  t,  where  x=x'+r<',  y—y',  z=z',  t=t', 
the  equation  satisfied  by  the  potential  is 


1/d,  c)\’_ 


d'<j> 

dx^'^dy-'^dz^ 


(4) 


or 


dxdt.  \c- 


1  \ 


This  equation  is  satisfied  by  the  potential  of  sources  of  sound 
in  motion  through  tlie  medium  with  uniform  velocity  v  in 
the  negative  x-direction.  It  is  also  the  equation  satisfied  by 
the  disturbance  velocity  potential  for  a  fixed  body  creating 
small  perturbations  from  an  oncoming  main  stream  of 
velocity  v  in  the  x-direction.  A  brief  derivation  from  liydro- 
(Kmamical  principles  is  given  in  appendix  A. 

It  is  known  from  tlic  classical  study  of  tlic  wave  equation 
(reference  16)  and  can  be  verified  by  direct  substitution  that 
a  solution  of  equation  (1)  is  transformed  to  a  solution  of 
equation  (4)  by  means  of  the  following  substitutions,  corre¬ 
sponding  to  a  combination  of  the  Lorentz  transfonnation 
and  a  Galilean  transformation: 


^  yT-dP 

y'=v 

z'  =  z 


xM 


vr=Ti^2 


(5) 


where  M,  the  Mach  number  of  the  main  flow,  is  »/c. 

For  the  purpose  of  studying  the  supersonic  case  (M;>l), 
it  is  more  cotivenient  to  employ  modifications  of  transforma¬ 


tion  (5)  obtained  by  multiplying  the  right-hand  side  by  the 
constant  l/'-y/l— or 


l-d/^ 

y 

■yl-dP 


V — t  -|- 


xM 


c{].-AP)  J 


(5a) 


The  particular  solution  of  equation  (4)  that  eorrespoiuls 
to  a  moving  source  will  be  seiui  in  the  following  discussion 
to  he  analogous  to  a  solution  of  equation  (1)  given  by  the 
sum  of  potentials  in  equations  (2)  and  (.3),  namely,  to 


The  desired  solution  of  ('quation  (4)  eon'esponding  to  equa¬ 
tion  (6)  is  obtained  with  the  aid  of  the  substitutions  (5a)  as 


4>o 

where 


-  [t  (‘  -  T  arli  -  0  M  '-iM  0]  (■> 


(The  term  in  ('filiations  (5a)  causes  no  difficulty 

since  only  the  s(|uari‘s  of  the  space  coordinate's  are  needed.) 
This  solution  for  may  ho  ('xpresse'd  in  the  form 

<t>,=  lU{t-r,)  +  J{t-T,)]  (7a) 

where 

_Al  X-?  r 
c'  A/^'-l+c 


_A1  X  — ^  _r 
c 

and  where  r  is  defined  as  in  equation  (8).  The  constant 
•‘Ul.’l.f)  could  of  course  have  been  includi'd  in  the  functional 
symbol  /  but  has  been  separated  for  convenience.  It  may 
lie  considt'i'ed  to  represent  the  space  variation  of  the  source 
strength  as  distinguished  from  the  time  variation  of  strength. 
For  a  moving  source  of  constant  strength  the  time  function 
may  be  considered  equal  to  unity  and  the  potential  ( xpressed 
us  (reference  2): 


It  will  be  recognized  that  the  solution,  equation  (7a),  is 
valid  in  a  conical  region,  the  so-called  “Mach  cone,’’  opening 
aft  of  the  moving  source.  Outside  of  this  conical  ri'gion, 
defined  by  the  equation  r=0,  the  flow  is  undisturbed. 
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The  result  expressed  by  equation  (7)  may  be  considered 
physically  from  two  points  of  view.  In  one,  as  considered 
by  Prandtl  (reference  2),  a  source  of  variable  strength 
moving  along  a  certain  path  is  replaced  by  a  continuous 
succession  of  fixed  source  pulses  distributed  along  its  path 
acting  consecutively  one  after  the  other.  Each  pulse,  con¬ 
sidered  fixed  in  an  absolute  coordinate  system,  emits  a 
spherical  wave  traveling  at  sound  speed  and  the  coordinates 
of  the  center  of  the  spherical  surface  are  ij,  f.  The 

radius  vector  R  of  a  point  (a;,  y,  z)  with  respect  to  this 
center  is 

B=  ^[x—  (y— rj)^+  {z—^Y 

The  time  at  which  the  spherical  wave  passes  the  point  {x,  y,  z) 
is 


c 


Eliminating  R  between  the  preceding  two  relations  results  in 

cH'^—(x—^—vtY—(y~riY—  (s— f)^=0 

The  roots  of  this  quadratic  equation  in  t  are  precisely  the 
quantities  ti  and  t2  defined  in  equation  (7a);  that  is,  the 
field  point  {x,  y,  z)  is  influenced  at  time  t  by  two  waves  which 
originated  at  times  r2  and  ti  earlier.  It  is  of  interest  to 
observe  that,  in  the  supersonic  flow,  both  roots  are  real  and 
positive  and  have  physical  significance;  whereas,  in  the 
subsonic  flow,  only  one  root  is  positive  and  of  physical  sig¬ 
nificance.  In  the  supersonic  case  the  field  of  influence  of  a 
source  is  the  particular  Mach  cone  with  vertex  at  the  souree, 
and  through  each  point  in  this  region  at  instant  t,  there  pass 
two  spherical  surfaces  representing  the  waves  originating  at 
times  Ti  and  t2  earlier  (fig.  2). 


Figure  2.— Field  of  influence  of  .■spherical  source  moving  at  a  constant  supersonic  velocity. 
Point  (X,  y,  2)  is  affected  at  time  t  by  two  pulses  originating  at  times  n  and  t2  earlier. 


I'rom  the  other  point  of  view  of  the  result  (equation  (7)), 
a  single  diverging  spherical  wave-pulse  is  considered.  Let 
this  wave  originate  at  the  point  (L  jj,  I)  at  a  time  T  (fig.  3) 
and  consider  its  effect  at  a  point  {x,  y,  z)  (within  the  Mach 
cone  whose  vortex  is  at  point  (L  17,  C))  moving  with  a  velocity 
greater  than  that  of  sound.  Clearly  at  a  later  time  T+ti 
the  moving  point  penetrates  the  wave  front  and  at  a  still 
later  time  T-l-r2  it  emerges  from  the  wave  front.  The 
potential  at  point  {x,  y,  z)  changes  only  on  entering  and  on 
leaving  the  wave  front  and  the  two  terms  in  equation  (7) 


Figure  3.— Influence  of  single  spherical  wave  pulse.  Spherical  wave  originating  at  points 
(f,  V,  r)  at  time  T  influences  point  (j,  y,  z)  fixed  relative  to  (^,  Vt  f)  and  moving  at  a  con¬ 
stant  supersonic  speed,  at  times  ti  and  n  later. 

correspond  to  these  two  effects.  The  factor  2  appearing  in 
the  potential  for  a  constant  source  moving  at  a  supersonic 
speed  also  has  its  origin  in  this  physical  fact,  in  contrast  to 
that  for  a  source  moving  at  a  subsonic  speed,  where  the  field 
point  penetrates  the  wave  front  but  never  emerges  and  where 
the  corresponding  factor  is  unity.  The  two-dimensional 
supersonic  case  involves  cylindrical  waves  and  the  potential 
of  the  point  (x,  y)  is  continuously  changing  from  the 
time  the  point  enters  to  the  time  it  emerges  from  tlic 
wave  (reference  13).  Observe  the  interesting  geometrical 
property  of  r  (equation  (8));  namely  2r  is  the  difference  of 
the  radii  of  the  spherical  wave  at  time  t2  and  at  time  ti,  that 

is,  r=-  (t2— Ti).  (Observe  also  that  the  potential  which  for- 

mally  appears  in  equation  (6)  as  the  sum  of  potentials,  half- 
advanced  and  half-retarded,  transforms  in  the  moving 
coordinates  to  a  sum  of  retarded  potentials  in  which  the 
original  retarded  part  is  associated  with  the  diverging 
spherical  concave  wave  from  which  the  point  is  emerging 
and  the  original  advanced  part  is  associated  with  a  diverging 
convex  wave  into  which  the  point  is  penetrating.)  Recent 
papers  of  interest  in  connection  with  moving  acoustical 
sources  are  references  19  and  20. 

SURFACE  DISTRIBUTION  OF  SOURCES 

Sources  and  sinks  of  the  type  will  now  be  distributed  to 
represent  the  upper  and  lower  surfaces  of  a  thin  airfoil. 
The  procedure  to  be  followed  is  that  used  in  the  two- 
dimensional  case  (references  10  to  13)  where  the  upper  and 
lower  surfaces  are  considered  separately.  Also  the  total  effect 
may  be  separated  into  an  effect  of  the  mean-camber  surface 
and  an  additive  effect  due  to  thickness  alone.  In  most  of 
the  applications,  unless  stated  to  the  contrary,  the  mean- 
camber  surface  is  considered. 

Let  a  continuous  distribution  of  sources  be  given  over  the 
mean-camber  surface.  The  airfoil  is  considered  so  thin  and 
flat  that  the  source  distribution  may  be  treated  in  the 
a;z-plane  (fig.  1).  The  airfoil  surface  may  be  considered 
moving  at  a  constant  speed  v  in  the  negative  x-direction  (or 
fixed  in  a  stream  moving  in  the  x-direction) .  The  effect  at 
a  point  (x,  y,  z)  at  time  f  of  a  distilbution  of  sources  of  posi¬ 
tion  magnitude  A{^,  0,  f)  is  given  by  an  appropriate 
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O 


iiil ('pillion  over  a  repoii  of  llic  ^f-plaiie  of  the  form 

4>{£,y,z,t)=f  f<t‘ad^  di  (9) 

w  ill'll'  <i>„  represi'iits  the  fiiuctioii  given  in  equation  (7)  with 

i;  =  0. 

Till'  total  ellVet  at  the  point  (x.  y,  z)  is  the  sum  of  the 
effeets  of  all  disturbances  having  their  origin  within  the 
Maeh  cone  with  vertex  at  point  (r,  ?/,  r)  and  opening  in  the 
upstream  direction.  This  conical  region  need  not  extend 
into  the  undisturbed  part  of  the  flow;  that  is,  it  need  not 
extend  beyond  the  most  forward  surface  envelope  of  the 
Mach  cones  of  influence  of  the  body.  There  are  essentially 
two  types  of  boundary  conditions  that  need  to  he  dis¬ 
tinguished,  designated  by  the  terms  “purely  supersonic.” 
and  “nii.xed  supersonic.”  A  point  of  the  hotmdary  belongs 
to  a  purely  supersonic  case  if  the  upstream  facing  Mach  cone 
contains,  in  the  part  of  the  rs-plane  not  considered  occupied 
by  the  body,  no  disturbed  fluid  having  a  component  normal 
to  this  surface.  Otherwise,  the  point  belongs  to  the  mixed 
supersonic  case,  A  sufficient  (but  not  necessary)  criterion 
for  the  purely  supersonic  case  is  that  the  component  of  the 
main  stream  normal  to  any  edge  or  contour  of  the  plan  form 
in  the  r^-planc  (contained  within  the  upstream  facing  Mach 
cone  of  the  given  point)  shall  be  stipersonic.  There  is  no 
downwash  ahead  of  the  body,  no  holes  arc  itt  the  body,  no 
spilling  of  fluid  occurs  around  edges,  and  no  reflecting  stirfaces 
are  in  the  flow  field.  In  this  case  the  upper  and  lower  sur¬ 
faces  of  the  airfoil  are  considered  to  act  independently  of 
each  other;  a  disttirbance  created  on  one  side  does  not  affect 
the  opposite  side.  The  flow  can  bo  considered  to  arise  from 
the  appropriate  movement  of  small  pistons  acting  at  the 
regulating  or  generating  surface.  This  condition  is  in  con¬ 
trast  to  that  of  the  mixed  supersonic  case,  for  which  the 
effect  of  the  disttirbance  spills  over  the  edges  or  sides,  and 
a  disturbed  fluid  region  (downwash)  may  exist  ahead  of  the 
body.  Thus,  points  of  a  triangular  surface,  moving  verte.x 
foremost  and  completely  outside  of  the  Mach  cone  associated 
with  the  vertex,  belong  to  the  purcl.v  supersonic  case.  If 
the  triangular  surface  is  inside  the  Mach  cone  associated 
with  the  vertex,  the  points  belong  to  the  mi.xed  supersonic 
case.  Of  course,  for  a  given  surface,  both  cases  may  be 
involved.  A  few  examples  are  shown  in  figure  4. 

In  the  purely  supersonic  case  the  circulation  concept  plays 
no  particular  role  and  the  drag  associated  with  lift  or  thick¬ 
ness  may  properly  be  denoted  as  wave  drag.  In  the  mi.xed 
case  the  flow  retains  subsonic  features  and  the  drag  associated 
with  the  lift  is  sometimes  denoted  as  induced  drag. 

.\lthough  the  treatment  given  for  the  purely  supersonic 
case  is  believed  exact  within  the  limitations  of  the  linearized 
tlieoiy,  an  exact  treatment  of  the  mixed  supersonic  case  is 
not  available.  These  problems  involve  greater  difficulties 
in  the  boundary  conditions,  for  the  flow  to  a  certain  extent 
acquires  features  of  a  subsonic  flow  in  that  the  fluid  field 
“senses”  the  approach  of  the  body.  Thus,  in  certain  cases, 
conditions  at  the  leading  edge,  at  the  trailing  edge,  and  in  the 
wake  must  be  specially  taken  into  account.  For  some  pur¬ 
poses  and  in  certain  problems,  however,  it  may  be  useful  to 
treat  the  mixed  supersonic  case  in  the  same  manner  as  the 


Figure  4.— Illustrations  of  plan  forms  with  “mixed  supersonic”  regions  (unshaded)  and 
“purely  supersonic”  regions  (shaded).  Mach  linos  are  shown  dashed. 


purely  supersonic  case  and  to  introduce  appropriate  correc¬ 
tion  factors. 

The  region  of  integration  in  equation  (9)  is  the  part  of  the 
body  (in  the  ^f-plane)  cut  out  by  the  upstream  opening 
■Mach  cone  with  vortex  at  point  (x,  y,  g).  This  region  in 
general  depends  on  the  plan  form  of  the  body  as  well  as  on 
the  point  (x,  y,  g).  With  the  understanding  that  the  leading 
point  of  the  body  is  at  i=0,  the  integration  may  bo  written 

'Pod^di  (10) 

where 

=  ^  ("o 

f2  =  gd-fo 

ii=x—yyM^-i 

The  limits  of  integration  f,  and  fo  in  equation  (10)  may  bo 
recognized  as  the  distances  from  the.  f-axis  to  the  near  and 
far  sides,  respectively,  of  the  hyperbola  defined  by  the  inter¬ 
section  of  the  cone  /•=0  and  the  plane  7j=0.  Thus,  from 
equation  (8),  with  7)  =  0,  f,  and  ("2  are  recognized  as  the  roots 
of  the  equation 

a/(F^T0GW)=O  (11) 

The  limit  in  equation  (10)  represents  the  f-coordinate  of 
the  vertex  of  the  hyperbola  and  is  defined  by  the  condition 
f,  =  f2,  that  is,  by  fo=0.  The  point  (^,,  z)  is  the  farthermost 
downstream  point  which  can  affect  the  point  (x,  y,  z). 
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BOUNDARY  CONDITION 

The  strength  of  the  distrilnitiou  of  singularities  in  equa¬ 
tion  (10)  will  now  he  determined  by  the  boundary  condition 
of  tangential  flow  along  the  airfoil  surface.  The  boundary 
condition  may  be  e.xpressed  as 


=  w(x,  z,t) 


bx^dt 


(12) 


where  the  airfoil  shape  is  defined  by  y—(j{r,  z,  t)  and  where 
the  two  terms  represent  the  normal  velocity  induced  by  the 
airfoil  shape  and  by  its  own  proper  motion.  It  is  shown  in 
appendi.x  B  and  can  also  be  made  clear  by  physical  reasoning 
that  as  y  approaches  zero  from  the  positive  side  (y-*+0) 

or,  briefly, 

A(,x,z)J{t)  =  -^;^^^j^^^wix,z,i)  (13) 

As  y  approaches  zero  in  the  negative  half  plane,  an  etpial 
and  opposite  result  is  obtained.  Equal  soiu-ce  distributions 
on  the  upper  and  lower  surfaces  therefore  result  in  a  discon¬ 
tinuous  vertical-velocity  distribution  near  the  plane  i/=0 
and  may  be  used  to  represent  symmetrical  thickness  distri¬ 
butions.  The  source  distribution  representing  a  thin  body 
with  arbitrary  thickness  distribution  is  in  general  unequal 
on  the  two  surfaces.  The  elFect  of  thickness  is  discussed  in 
a  separate  section.  A  representation  of  the  mean-cambiu' 
surface  alone  may  be  obtained  by  placing  equal  and  opposite 
sources  on  the  under  surface  in  pro.ximity  to  the  sources  on 
the  upper  surface.  The  potential  <t)  is  to  be  understood  in 
the  subsequent  analysis  to  be  prefi.xcd  by  a  ±  sign,  plus  for 
the  upper  surface  and  minus  for  the  lower  surface.  The 
vertical  velocity  will  in  general  be  measured  positive  upward. 

It  is  convenient  to  express  the  vertical  velocity  in  equa¬ 
tion  (13)  in  separated  form 

w{x,z,t)  =  \V{x,z)w{t)  (14) 

where 

W(x,z)  =  —  27r(ilP—  l)A(x,z) 
w(,t)  =J{t) 

SURFACE  POTENTIAL 

The  total  potential  for  ?/=0  may  now  be  expressed  by 
means  of  equations  (10)  and  (14)  as 

4>(x,z,t)  =  j~j^  4>o  df 


W(U) 


W(<  — T,)-l-W(<  — To) 


dr  d^  (15) 


where,  for  y=0  (sec  equations  (7a),  (10),  and  (11)), 

=:  tII  _  \'(r— ri)(h— r) 

cl3-  c/3 

^!(x—i) .  A (r — ri)(r2‘~r) 

ri=3 — To 

S'2=3+S'o 

;3=  vil/'-T 

and  w'here  it  is  understood  that  \V(^,  r)=d  at  any  pi 
the  body  or  where  the  integrand  is  not  real. 

Equation  (15)  may  be  put  into  a  simpler  form  by 
tution  of  a  new  variable  B  instead  of  r.  which  is  ob 
from  the  relation  (see  appendix  B) 

2r=(r2~ri)  s-i-rj+ri 

or 

r=ro  cos  0+x 

The  surface  potential  (equation  (15))  may  then  be  writ 

<t>{x,  z,  0  =  -2“-jsj]  r{9)][w(f-T,)+w(f-T2)]  I 

where 

ri=^2-  (Af-sin  B) 
t2=^^y  (M-l-sin  B) 

0=cos"‘  ^ — ?  0 

Equation  (15)  represents  the  central  result  of  the  ai 
and  within  the  limitations  already  discussed  may  bo  a 
to  wings  of  any  plan  form  in  steady  motion  or  perfo 
small  oscillations.  In  the  stationary  or  steady  case,  < 
not  depend  on  time  and  the  function  wit)  is  to  he  rc| 
by  unity.  Then,  in  equation  (15),  w{t—Ti)+w{t—T2: 
be  replaced  by  2. 

PRESSURE  RELATIONS 

For  the  sake  of  reference,  relations  for  the  pressui 
the  lift  and  drag  forces  are  given  here.  The  distur 
pressure  (local  static  pressure  minus  the  pressure  i 
undisturbed  stream)  may  be  written  as 
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'Plio  pi'ossm'c  (lifl'eronco  (posilivc  if  acting  ilownward)  at  any 
])oint  (j:,  r)  may  l)c  expressed  as 

Ap  =  i>v—l>L 

wliere  the  subscripts  l~  and  L  refer  to  the  upper  and  lower 
surfaces.  For  the  mean-camber  surface  =  uml 

A _ O-  t  /.TA 


'Phe  total  forces  on  the  airfoil  in  the  (/-direction  and 
a-dircction  arc  given  by 

F=Lift=y’y’2i  d£  dz 

A'=Drag=  —ffp  dy  dz 

where  the  integration  is  to  be  taken  over  the  eomplete  airfoil 
surface.  Expressed  as  integrations  over  the  plan  form 

F=  JJ  (ii^—iic)dxdz 

It  is  often  convenient  to  sejiarate  the  slope  terms  as  follows: 


The  integration  with  respect  to  6  may  be  readily  performed 
with  the  aid  of  the  relation 


where  a  is  the  conventional  direction  of  the  main  stream  with 
respect  to  a  reference  chord,  and  <rc;  and  ar.  arc  I  he  local  slopes 
of  the  airfoil  surfaces  measured  with  respect  to  the  reference 
chord  and  positive  in  the  same  sense  as  a. 

APPLICATIONS 

WING  OF  INFINITE  SPAN  AND  ZERO  SWEEP 

For  the  first  application  of  equation  (1,5)  the  results  for 
both  the  oscillating  and  steady  two-dimensional  case  will  be 
derived.  For  the  harmonically  oscillating  wing  having 
identical  motion  in  every  chordwisc  section,  the  vertical 
velocity  can  be  written  in  the  complex  form 

■w{x,  t)  =  ir(x)e‘‘“‘ 

Then 

IP  (f  —  r  1 ) -p  w  ( t  —  Tj)  =  e ‘  ~ 


Equation  (15a)  becomes 


2e  ~  cos 


4,(x,  t)  =  J^'fos  ^  sin  e)  d9di 

_  (19) 

where 


Finally 


—  (  cos  (X  sin  d)dd  =  jQ(},) 

TT  Jo 

fiiut  f*X 

0(x,f)  =  -y  W(i)lii,  x)d^ 


This  result  for  the  velocity  potential  is  identical  with  equa¬ 
tion  (11)  of  reference  13  and  is  used  therein  as  a  basis  for 
calculation  of  the  nonstationary  two-dimensional  case. 

In  the  steady  case,  a)  =  0  and  /(J,  x)=l.  The  expression 
for  the  velocity  potential  is 

0(x)  =  -ij^r(^)d«  (22) 

(iv 

where  This  formula  or  the  pressure  relation 

dif)  dy 
v=  —  ou  .  =0  -T-  -r 


applied  to  both  the  upper  and  lower  surfaces  of  the  airfoil 
leads  to  all  the  results  of  the  Ackeret  theory. 

WING  OF  INFINITE  SPAN  WITH  ANGLE  OF  SWEEP 

Consider  an  infinite  wing  with  angle  of  sweep  A  (fig.  5),  and 
assume  that  all  sections  in  the  flight  direction  arc  identical 
in  shape  and  that  the  wing  is  undergoing  harmonic  motion. 
In  general,  the  vertical  velocity  w  can  bo  written  in  the  com¬ 
plex  form 

w(x,  2,  t)  —  W{x,  2)e'"‘ 


Figure  5.— Sketch  for  wing  of  infinite  span  wit  h  angle  of  sweep  showing  region  of 
integration  (shaded). 
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]f  each  section  normal  to  the  leatlinf;  edfje  is  performins;  the 
same  motion,  the  form  of  Thf-'',  s)  is  W(x—2  tan  A).  If  the 
wins;  is  assumed  to  perform  pure  vertical  motion  alone,  then 
2)  is  a  constant.  If  the  wing;  is  assumed  to  rotate  about 
an  a.xis  j'=Constant,  then  IFf-r,  s)  is  of  the  form  HXt). 

The  potential  is  of  the  form  (fig.  5) 

I'l  r-fr  re, 

ij>{x,z,t)=  F(i6(l^—  I  F  (19  di  (23) 

J(i  Jn  Jjj  Jo 


F(i,  6, 0  =  -0^  f  (®))  + 


and  where 


1— (3  cot  A 


1  +  /3  cot  A 


fli  =  cos'‘ 


cot  A— 2 


The  values  of  the  limits  ^3  and  are  found  by  solving  for  f 
in  the  relations  s"i  =  |  cot  A  and  cot  A  which  represent 
the  intersections  of  the  Mach  lines  through  x  with  the  lead¬ 
ing  edge.  The  limit  6=6i  corresponds  to  cot  A,  the 
leading-edge  line. 

When  ir(|,  2)  is  a  constant  or  a  function  of  (  only,  the 
velocity  potential  can  be  c.xpressod  as 

0(*,  2,  <)=-y  x)d«- Jj*‘ir(«/,({.  X,  2)d?] 

(24) 

whore  /(J,  x)  is  as  defined  previously  and 


1  ,  x-(  .\r  re, 

/.(^,x,2)=^e-'“-^  ( 


Observe  tliat  the  integral  involved  in  equation  (25)  for  /, 
reduces  to  the  Bessel  function  of  zero  order  when  $,  =  ir  as 
in  equation  (19).  This  interesting  integral  may  therefore 
be  called  an  “incomplete”  Bessel  function  of  zero  order. 
Systematic  investigation  of  its  properties  would  appear  to 
he  desirable. 

For  the  infinite  swept  wing  in  the  steady  case  the  frequency 
(j>  may  he  made  equal  to  zero  in  equation  (23).  Consider 
as  a  simple  e.xamplo  the  case  of  a  thin  wing  at  a  small  constant 

angle  of  attack  a.  that  is,  — «•  het  the  angle  of  sweep 

be  less  than  the  complement  of  the  Mach  angle,  that  is, 
(3  cot  A>1.  (Otherwise  the  case  involves  the  mi.xed- 
supersonie-llow  eondilions.)  li'rom  equation  (24)  with  a)=0, 

/(J,  x)  =  l,  and  /,(J,  X,  2)  =  ®‘-, 

TT 


The  local  pressure  difi'erence  is  given  by 

^_2pr*a  0  cot  A _ 

/3  V/3- coC  A-1 

This  equation  reduces,  for  A  =  0,  to  the,  Ackerel  resu 


Let  the  inde.x  n  refer  to  quantities  measured  norma 
leading  edge.  Then 

a„  =  a  sec  A 
v„=v  cos  A 

’‘=M  cos  A 
c 

and 

2pv„^an 

a  result  similar  in  form  to  the  e.xpression  for  po  and 
slated  by  Busemann  (reference  17)  in  1935.  (S. 

reference  61) 

The  harmonically  oscillating  case  with  ir(x,  z)  a- 
t  o  be  of  the  form  \V  (x—z  tan  A)  leads  in  a  similar  mas 
a  result  analogous  to  equation  (20). 

RECTANGULAR  WING  OF  FINITE  SPAN  (ZERO  SWEEP) 

Consider  a  harmonically  oscillating  rectangular  \ 
finite  span  as  in  figure  6.  Region  I  is  described  ns 
supersonic  and  region  II  as  mi.xed  supersonic.  The 
the  aspect  ratio  and  the  stream  Mach  number,  the  rel 
smaller  the  region  II  becomes. 

The  potential  for  region  I  for  identical  motion  ( 
chordwise  section  is  e.xnctly  that  given  for  the  infinii 
in  equation  (20) ;  however,  more  general  types  of 
involving  spanwise  variation  may  also  be  treated 
example,  let  the  wing  perform  harmonic  oscillati' 
vertical  bending  and  in  torsion  about  a  spanwise  a.xi 
in  certain  prescribed  spanwise  modes.  Then,  with 
h  used  to  describe  angle  of  attack  and  vertical  p 
(fig.  7), 

a  =  oci{z)a2(t) ) 

h=hi(z)h2(t) ) 

where  a,(z}  and  h,(z)  represent  spanwise  modes  and 

a2(t)  =  ane‘“‘ 
h2{t)=hoe'“‘ 

and  ao  and  ha  are  constant  complex  amplitudes.  The  v 
velocity  (w  measured  positive  upward,  h  positive 
ward)  may  be  expressed  as 

w{x,  z,  t)  =  —  [ra-|-fi-l-(x— x„)a] 
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> 


FiGVRE  7.— Sketch  tor  oscilhitintr  rcctuiigular  wing  (A  Is  positive  downward, 
a  is  positive  clockwise). 


Lot  (ho  potoulial  (oquation  (15))  ho  soparatod  into  (lio 
form 


(30) 


where'  the  various  0’s  are  associated  witli  the  correspond intr 
variables  in  equation  (29). 


With  the  use  of  equations  (15a),  these  potentials  may  bo 
expressed  as 

J  Q:,(r)  <^os  sin  9^  1 

fi,  cos(^^^^sin  >  (31) 

0„=“^  (^— ri,)Q:,(i-)  cos^-^^sin 

where 

_ uAf  (x — J) 

and,  expressed  as  a  function  of  6, 

“i(f)  =  “i(3+fo  eos  B) 

7t,(f)=/q(3  +  fo  cos  9) 

If  the  modal  functions  in  equations  (31)  are  a,=/!i  =  l,  the 
potential  corresponds  to  that  given  by  equation  (20)  for  the 
two-dimensional  case.  (See  also  equation  (14)  of  reference  13.) 
It  is  of  interest  to  consider  modal  functions  for  ai  and 
hi  of  the  type  (f/.f)’’  where  s  is  the  semispan.  For  modal 
functions  of  this  form  the  typical  integral  involved  in 
equations  (31)  may  bo  expressed  as 

F,= (s-t-fo  cos  9)"  cos  sin  b'J  dB  (32) 
With  the  substitution  of  ^  —  9  for  B,  F„  may  be  written  as 


F„=  (2-i- fo  sin  B) "  cos  cos  9^  dB+ 

(s— fo  sin  9)"  cos  (^-^cos  9^  dB 

The  further  reduction  of  F„  is  made  with  the  aid  of  the  fol¬ 
lowing  relation  (reference  21) : 

TV  cos  (X  cos  9)  sin“  9  dB 

2/ 

For  example,  the  case  n=0  corresponds  to  constant  modes 
and  yields  for  the  potential  the  result  already  given  by 
equation  (20).  The  case  ri.=  l  corresponds  to  linear  modes 
and  the  function  Fi  becomes 

This  relation  utilized  in  the  equation  for  the  potential  yields 
a  result  that  is  the  two-dimensional-case  result  multiplied 
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by  the  factor  z/s.  The  case  n=2  corresponds  to  parabolic 
modes  and  the  function  F,  becomes 

IMien  Fi  is  used  in  equations  (Rl),  the  Jn  term  yields  an 
integral  of  the  type  given  by  equation  (20).  With  the  use 
of  the  relation  Ji(X)  =  —  Jo'(^),  the  Ji  term  also  yields  an 
integral  of  the  same  form.  This  type  of  reduction  to  the 
form  of  equation  (20)  may  bo  made  in  general  for  any  integral 
index  n  by  means  of  the  recurrence  formulas  for  Bessel  func¬ 
tions,  and  thus  use  may  be  made  of  the  numerical  procedures 
used  for  equation  (RO).  (See  reference  IR.) 

It  may  be  of  interest  to  treat  the  potential  for  the  mi.xed 
supersonic  region  II  (lig.  6)  as  though  it  were  part  of  a  purely 
supersonic  region.  The  equations  corresponding  to  equa¬ 
tions  (31)  are 

Jo"  "'‘'Jo  Cv/ 

(m 

and  similar  equations  for  0/i  and  (/>i.  The  limit  is  found 
as  the  value  of  i  for  which  f2=s  or 

^j=x-(s— z)j3 

The  limit  9=0  corresponds  to  r=f2  and  the  limit  9=92 
corresponds  to  !■=«  or,  from  equation  (15a), 

9,=cos- 

The  last  term  in  equation  (RR)  leads  to  integrals  of  the 
“incomplete”  Bessel  function  type  as  mentioned  for  the  ease 
of  the  infinite  wing  with  angle  of  sweep. 

The  foregoing  results  for  the  oscillating  rectangular  wing 
will  now  be  specialized  to  the  steady  case  (a)=0,  2=0, 
oi(r)=ii  “2(0  =  a,  the  constant  angle  of  attack).  Then, 
from  equations  (31),  the  velocity  potential  for  region  I,  is 

X  (34) 


It  may  be  observed  that  y  is  constant  along  rays  fi' 
tip  Constant.  Along  the  ray  corresponding 

5  —  z 

Mach  line  from  the  tip,  /3=  1  and  p  takes  on  the  C( 

value  Pi.  iVlong  the  ray  corresponding  to  the  tip  £= 
of  this  value  is  obtained.  This  edge  condition  is  pliA 
incorrect  since  the  assumption  of  the  independence 
two  surfaces  of  the  airfoil  is  not  correct  near  the  tip. 

This  particular  problem  has  been  treated  by  Bin 
(reference  4)  by  his  method  of  conical  or  perspcctiv 
metry.  The  condition  along  the  ray  corresponding 
tip  is  p=0  and  Busemann's  result  for  region  II  is 


Tii= 


2pv'‘a  1 
d  jr 


The  total  lift  over  region  II  is  one-half  of  that  of  at 
area  of  region  I.  A  comparison  of  this  result  and 
tions  (36)  is  shown  in  figure  8.  This  comparison  gi 
indication  of  the  errors  involved  in  the  assumption  0 
pendcncc  of  the  two  surfaces  near  the  rectangular  t 
conversely,  it  gives  an  indication  of  the  appropriate  ■ 
tion  factors  required  to  allow  for  the  tip  effect.  It  a 
that  equatio'ns  (36)  overestimate  the  lift  over  all  of  rci 
2 

by  a  factor  1—-  or  by  appro.ximately  36  percent.  Ac 

equations  (36)  apply  to  the  edge  of  a  rectangulai 
adjacent  to  a  straight  surface  barrier  at  zero  angle  of  1 

THICKNESS  DISTRIBUTION 

It  has  already  been  remarked  that  the  treatmei 
ployed  in  the  analysis  mainly  for  the  mean-camber  s 
can  also  be  applied  to  obtain  the  effect  of  thicknes 
equation  (15)  the  vertical  velocity  IF(f,  f)  may  bo  sp 
for  both  the  upper  and  the  lower  surfare. 

jVs  an  example,  consider  a  plan  form  such  as  that  shi 
figure  9  in  steady  supersonic  flow.  Let  the  airfoil  s 
shape,  for  convenience  chosen  symmetrical  and  indepi 
of  span,  be  defined  in  the  center  section  by  y=g{x)  (•■ 
any  other  section  by  y=tj(x—z  tan  A)).  Then,  for  the 
surface, 

IF(x,'^^7)=ea-f 


For  region  II,  from  equation  (33), 


and,  for  the  lower  surface. 


va  Va  r^’  /s  —  Z  -A  /..-N 

The  actual  integration  in  equation  (35)  may  he  easily  per¬ 
formed  but  is  not  required  for  the  purpose  of  obtaining  the 
local  pressure. 

The  local-pressure  difference  is  directly  obtained  for 
regions  I  and  II  from  equations  (34)  and  (35)  as 


lT(x.  z)=pa-i-g' 

where  (/  is  the  derivative  of  g  with  respect  to  its  argumi 
The  velocity  potentials  in  the  various  regions  in  fq 
are  of  the  form 

nr  r(i  rs,  N 
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Figure  I’ressure  for  rectangular  edge  as  given  by  Husemami  and  by  llic  approximation 
given  io  equations  (30). 


Figure  y,— Skuicb  for  swept  wing  for  the  case  /1>1  showing  regions  I,  II,  III. 


The  limits  in  the  foregoing  integrals  are  as  follows:  the 
limits  and  6=9^  correspond,  respectively,  to  the  leading- 
edge  linos  1"  =  ^  cot  A  and  f=  — ?  cot  A:  9=0  and  9=x 
correspond,  respectiv'ely,  to  the  Mach  lines  f=s'2  and 
1  =  Si;  t3  m**!  fs  Ric  obtained,  respectively,  from  the  relations 
fi  =  J  cot  A  and  u  =  —  ^  cot  A;  9=9-2  corresponds  to  f=s;  is 
obtained  from  the  relation  s'2=«;  it  is  obtained  from 


=  {  cot  A;  and  is  the  value  of  J  for  the  leading  edge  of 
the  tip.  Then 


—  t  cot  A— g 

■c  — ? 


) 


92  =  cos  ‘  jS 

X  s 

_x—z^ 

y  X-1-2/3 

z— g/3 
-‘“l-fM 

{2=z-(s-s)/3 

^5=«  tan  A 

A=^  cot  A 

If,  for  example,  the'distribution  function  F  is  a  constant  K, 


<I>,=K 


X  cot  A— 
A^-1 


z 


nja  Fj{, 

<f>iii  =  0i~~r  (ir  — 9i)rfJ+  r  (it  — Sjrff 
irjfi  Jo 


The  corresponding  local  pressures  are 


pvK  A 


poK  A  /,  1 

P  JaF=1\  a+b) 


Piii= 


pvK  A 


/3  V^^- 
pvK  A 


1  ir\ 


gQo-i  — — 


where 


■  V  i-c^ 


~  —  tan  ■-%/ 

•  1  TT  V 


A=j3  cot  A 


B= 


{s—z)0 

_ o  fori  A 


The  constant  K  may  be  interpreted  as  va  associated  with 
constant  angle  of  attack.  In  this  case,  region  II  is  to  be 
regarded  as  a  mixed  supersonic  region  and  the  result  given 
is  not  the  appropriate  solution  for  this  region.  If  the  con¬ 
stant  K  is  interpreted  as  £i'=Constant,  the  results  are 
applicable  to  a  thin  symmetrical  wedge  of  half  vertex  angle 
K  and  may  be  employed  to  yield  the  wave  drag  according  to 
the  linearized  treatment. 
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Jones  (reference  5)  treats  symmetrical  airfoils  of  various 
plan  forms  at  zero  lift  by  use  ot  pressure  potential.  The 
use  of  velocity  potential  leads  to  the  same  results  as  given 
in  referenee  o.  'riuis,  ecpiations  (V.i)  and  (14)  of  reference  5 
for  a  wedge  correspond  to  the  preceding  results.  'I'hc 
velocity  potential  in  general  is  more  useful  to  treat  pressure 
distributions  for  a  given  body;  whereas,  the  pressure  poten¬ 
tial  may  be  more  readily  adapted  to  treat  airfoil  shapes  and 
plan  forms  associated  with  desired  types  of  distributions 
of  pressure. 


steady  case  of  a  vanishingly  thin  surface  at  angle  of  ai 
tile  velocity  potentials  and  pressure  relations  for  it 
and  III  are  equivalent  to  those  just  discussed  in  the  pr 
section.  The  lifl,  AL  on  a  strip  Ax  of  the  triangle,  loc 
abseissa  x  from  tlu'  v(>rte.\  is  given  by 

/*j:  cot  A 

AL=Ax  Ap  ds 

J  —X  cot  A 


_r  ipv'^a 


L  ff- 


+- 


TRIANGIII.AII  PLAN  FOHM 


The  triangular  wing  (fig.  10)  e.xteuding  across  the  Mach 
lints  from  the  verte.x  may  serve  as  a  liual  e.xample.  For  the 


where  tlu'  two  tiTins  correspond  to  the  integral  ioi 
regions  I  and  III,  respectively,  and  where  A=^  cot  A. 


FiOURS  10.— I'riuiigiilar  wIiik  in  a  siupcrsunic  stream  and  pressure  <Ustrlbution.  ('use 
sketched  corresponds  to  A^30®,  A/<=V2. 


AL—  Ax 

The  area  of  the  strip  is  2x  Ax  cot  A,  and  hence  the  lift 
cient  is  independent  of  x  and  equal  to 

p  _ _ _ 

/I  M2AxAx\ 

J 

Gurevich  (referenee  9)  treats  this  case,  and  his  rclatii 
bo  shown  to  be  equivalent  to  the  foregoing  ones.  Th 
sure  distribution  is  illustrated  in  figure  10.  where 
reference  pressure,  is  2pv^a/^.  Observe  that  the  p 
area  above  the,  unit  ordinate  cancels  the  area  of  p 
deficiency  below  the  unit  ordinate.  Also  shown  in  fi;. 
is  the  distribution  of  pressure  as  the  half  vertex  angle 
triangle  approaches  the  Alach  angle. 

The  triangular  wing  inside  the  Mach  cone  from  the 
requires  a  more  elaborate  treatment  (references  7  to  ' 


Langley  Memorial  Aeronautical  Laboratory, 
iVational  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va.,  June  4,  1947. 


APPENDIX  A 

DIFFERENTIAL  EQUATION  FOR  THE  VELOCITY  POTENTIAL 

A  derivation  of  equation  (4)  is  given  briefly  here.  The 
condition  for  irrotational  flow  is 

curl  v=0  (Al) 

and  this  relation  implies  that  a  scalar  velocity  potential  <j> 
e.xisls,  such  that 

v  =  grad  <t>  (A2) 

The  general  etiuation  of  continuity 

^+divpv=0 

may  bo  written  as 

(A3) 

where  differentiation  following  the  particle  is  denoted  by 
D  d  ,  ,  , 


With  the  aid  of  these  two  relations  the  first  term  in  equa¬ 
tion  (A3)  becomes 


0 _ _  1 

(■-  \W 


For  small  perturbations  from  the  main  stream  of  velocity  r 
in  the  s-direction,  c  may  be  considered  equal  to  the  constant 
spoetl  of  sound  in  the  undisturbed  medium  and,  in  compari¬ 
son  with  V,  v,=0,  and  Vi=i'.  Then 


c-  \  dt^  dxdt^ 


"  dx-) 


With  this  relation  used  in  equation  (A3)  the  equation  for 
the  velocity  potential  may  be  put  in  the  form  giv'en  in  equa¬ 
tion  (4)  of  the  analysis. 

APPENDIX  B 

EVALUATION  OF  m 


and  V^=div  grad  is  the  Laplacian  operator. 

From  Euler’s  equations,  or  from  the  general  Bernoulli 
relation, 


where  a  space  constant  function  of  time  has  been  included 
in  <t>,  and  where  it  has  been  assumed  that  p  is  a  function  of 
p  only.  With  the  use  of  equation  (A4)  and  the  acoustic 
relation, 

,  dp 

where  c  is  the  local  variable  speed  of  sound,  it  follows  that 


=  --  grad  p 


In  order  to  determine  the  limit  of  ^  as  y->0,  it  is  conven¬ 
ient  to  make  use  of  the  following  substitution: 

2^=  (f'j— Zi)  cos  (Bl) 

The  expression  for  0  (equation  (10))  may  be  written  with  the 
aid  of  the  following  relations  (see  equation  (11)): 

-5T™* 

S'o  sin  0 

<f>=/3  £  ^l(?,0,2+fo  cos  6){Ji+Ji'}dedi  (B2) 


I 

i)t\U^2)  p 
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By  the  rule  for  differentiation  of  a  definite  integi-al, 
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1 

Pdy  dy 


dy 


rj> 


(B3) 


Make  use  of  the  following  relations: 

dYl_ _ 3.4 _ 3(2+ j'o  cos  0) 

dy  ~3(2-(-fo  (R)s>)  3i/ 

_3+3fo. 


^(A+J=)_ _ 

dy  3(«-t,) 


^  cos  8 
dz  dy 

a/,  3«-r,) 


+ 


df,  3(<-t,) 


dy  ^3(f— To)  dy 


Then,  by  integration  by  parts,  the  next  to  the  last  integral  in 
equation  (B3)  becomes  ^withcosddd=(/y,  ^  (/i+/2)=a^ 

[<■'''+-*>  If ’“’'"5  Ml 
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(/i  — /s)  sin  9  cos  9— 


A  ^ 

e/3^  dt‘ 


2  (/1+/2)  cos-9jc/9(f{ 


where  the  first  term  vanishes  because /i=/2  at  9=0  and  8=ir. 
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-2t(+P-  l)+(z,0,2)/(<)  (B4) 


Since  g-  changes  sign  as  y  changes  sign,  it  follows  that 

as  y  approaches  zero  from  the  negative  side  an  equal  and 
opposite  result  is  obtained. 
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Positive  directions  of  axes  and  angles  (forces  and  momenta)  are  shown  by  arrows 
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Pitch 
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X - >Y 
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Angle  of  set  of  control  surface  (relative  to  neutral 
position),  5.  (Indicate  surface  by  proper  subscript.) 


Absolute  coeflScieuts  of  moment 

n-  ^  r  -M. 

(roUmg)  (pitchmg) 


JV 
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(yawing) 


4.  PBOPELLER  SYMBOLS 


D  Diameter 

p  Geometric  pitch 

pjD  Pitch  ratio 

V'  Inflow  velocity 

V,  Slipstream  velocity 

T  Thrust,  absolute  coefficient  Ct= 

Q  Torque,  absolute  coefficient  Cq— 


!  hp=76.04  kg-m/s=550  ft-lb/sec 
1  metric  horsepower=0.9863  hp 
1  mph=0.4470  mps 
1  mps=2.2369  mph 
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.  ^ 
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pn^IP 

Power,  absolute  coefficient  Cp= 


Speed-power  coefficient  = 


7: 


pV^ 

Pn^ 


Elfficiency 

Eevolutions  per  second,  rps 

.  ...  _  .  _  /  V  \ 

Effective  helix  angie^an"'^g - j 


5.  NUMERICAL  RELATIONS 


1  lb= 0.4536  kg 
1  kg= 2.2046  lb 
1  mi=  1,609.35  m=5,280  ft 
1  m=3.2808  ft 
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